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ABSTRACT Images produced from the Maximum Correlation Meth- 
od sometimes suffer from visible striping artifacts, especially for areas of 
extended sources. Possible causes are different baseline levels and calibra- 
tion errors in the detectors. We incorporated these factors into the MCM 
algorithm, and tested the effects of different constraints on the output im- 
age. The result shows significant visual improvement over the standard 
MCM method. In some areas the new images show intelligible structures 
that are otherwise corrupted by striping artifacts, and the removal of 
these artifacts could enhance performance of object classification algo- 
rithms. The constraints were also tested on low surface brightness areas, 
arid were found to be effective in reducing the noise level. 


INTRODUCTION 

The 1PAC program LA UN DR invokes several one- dimensional flat fielding and 
deglitching techniques. For the purpose of dcstriping the one-dimensional algo- 
rithm works well for regions with a well-defined baseline, but the results are not 
satisfactory for regions where structure exists at all spatial frequencies (Fowler 
and Melnyk 1990). 

Another IPAC utility KESTER, developed by Do K ester, is similar in prin- 
ciple to the approach we take. The idea is to process the data with YORIC to 
a certain iteration to obtain an image, which is then used to simulate a set of 
detector flux measurements. The original data are then calibrated against the 
simulated ones. 

Our approach is to combine image construction and the dcstriping process. 
Since the striping gets amplified through the iterations, the idea of applying 
constraints to the correction factors is natural. We investigated the possibilities 
of both additive and multiplicative offsets, and from our preliminary analysis it 
seems the multiplicative treatment is superior. 

ALGORITHM 

The Maximum Correlation Method can be derived based on a Maximum Like- 
lihood scheme. Assuming a Poisson process, the likelihood function for a set of 
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measurements D{ is 
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where F t is the noiseless measurement value. If r,y is the value of the ith foot- 
print’s response function at image pixel /j, then 
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To maximize the likelihood estimate for /j, we have 
d In Q 

-JJ- = ^ + A ln f i - A ln A + A] = £1-^ + A/-Fi] = 0. (3) 

From this an iterative algorithm can be built (Aumann, Fowler and Melnyk 
1990), 
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starting from a uniform image = /£ = for all j. 

For /RAS data however, differing baselines and calibrations exist among the 
detectors, which results in striping artifacts in the HiRes images. Figure 1 (a) 
and (b) show the 1st and 20th iteration HiRes images for a field in the p Ophi- 
uchus region at 100 /xm. The input data were processed with the IPAC utility 
LAUNDR 1 , which applied one dimensional flat fielding to the flux measurements 
(Fowler and Melnyk, 1990). 

We now incorporate the estimation of leg offsets with image construction, 
i.e. we’ll first try to maximize the likelihood function by choosing proper offsets 
for the legs, then proceed to compute the correction factors from the modified 

ZV s. 


A, Assuming the offset is additive, 

D* = - Oi is the corrected measurement value. Oi is the offset, which is 

the same for all V s within the same scanline X. 
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5 The LAUNDlt’ing was carried out for a 6.5° x 6.5° field, which was then clipped to obtain this 
1° x 1° region. 
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Maximizing P with choice of Ol, 


d\nP 

0O L 


E [-fi+(A--Oi)lnCi 

VOL ; i„ leg i 

(Di - 0 L )\n(D t - 0 L ) + (D t - 0 L )} 

£ [-lnfi+ln(A-0 L )] 

i in leg L 

o 


(7) 


So we have, 

J] n C * = 1 - (8) 

i in leg L * i in leg L 

This polynomial equation in Oi can be solved with Newton’s method. Ol s 
are computed starting from the second iteration, since in the first iteration 
f9 - Vi j f° = const, and there is no reason to regard it as an approximation 
of the true i\. 


B. Assuming the offset is multiplicative. 

Due to calibration error in the gains and the quasi-logarithmic digitization used 
in IRAS, the offset could be multiplicative in nature. In this case D * = miD l 
and din P/ dm i = 0 gives 
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RESULTS AND DISCUSSIONS 

The test field we used is the 100 /an band of a region in /> Ophiuchus (RA 
161i23m, Dec -23°30 / ). It covers a 1.07° x 1.07° field (257 x 257 15 arcscc 
pixels), with a total of 12112 footprints in 227 scanlines. 

The 20th iteration image produced from method A is shown in Figure 2(a). 
The calculated leg offsets are roughly Gaussian. 

Corresponding image from method B is shown in Figure 2(b). The calcu- 
lated leg offsets have a standard deviation of 0.12. 

Both Figures 2(a) and 2(b) show significant improvement over Figure 1(b). 
Striping artifacts are greatly reduced, while believable high spatial frequency 
features are retained. Figures 2(a) and 2(b) have much smoother contours than 
Figure 1(b). 

We also notice that Figure 2(b) is smoother than Figure 2(a) in the high 
flux region. We assume the reason is that the additive offsets calculated were 
not able to compensate for the differences among nearby legs in the high flux 
region. 




FIGURE 2 20th iteration, additive and multiplicative offsets 
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FIGURE 3 Recovery of Artificially Introduced Offsets 


We carried out a simple statistical test oil the nature of the leg offsets by 
comparing the likelihood function P obtained from the two assumptions for the 
first calibration. The multiplicative assumption results in a bigger likelihood for 
145 legs out of a total of 227. This again seems to indicate the multiplicative 
assumption is closer to the truth. 

In addition, we tried the constrained algorithms on the 100 /im band of 
M101. This field features intersecting scanlines, and because of the presence of 
a well-defined background, is calibrated quite well with the LAUNDR program. 
However the constrained image shows a reduction in noise. (Figure not shown. 
Calculated offsets have a standard deviation of 0.05). 

VERIFICATION 

We tested the algorithms’ capability to recover offsets artificially introduced to a 
set of simulated detector fluxes. The detector fluxes were computed from a first 
iteration image of the same region, which is virtually stripe free. Gaussianly gen- 
erated offsets (additive or multiplicative, magnitude being equal to that found 
for actual measurements) were then applied to the legs. This made-up set of 
detector fluxes was fed to our program, and scatter plots of the recovered vs. 
introduced offsets are shown in Figure 3. (left: Y-axis: recovered additive offset, 
X-axis: introduced additive offset; right: Y-axis: log of recovered multiplicative 
offset, X-axis: log of introduced multiplicative offset.) 

In the additive case, we found the residual offsets after attempt of recovery 
have standard deviation roughly one-sixth of that for introduced offsets, meaning 
a 6 fold decrease in the striping amplitude (or 15 dB decrease in power). The 
performance is limited by the density of the scanlines. 

In the multiplicative case, the introduced offsets were a Gaussian random 
variable with standard deviation 0.12. The standard deviation of residual offset 
is 0.024. 

In these tests the recovered images are visually indistinguishable from the 
noiseless input image. This suggests that the algorithms are capable of doing 
what they are supposed to. So the residual striping seen in Figures 2(a) and (b) 
is due to inaccuracy of the underlying assumptions about the nature of the off- 
sets. Test runs on neighboring fields showed indeed the offsets vary considerably 
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within each leg over the spatial scale of 1°. One of the many possible causes 
is the change in responsivity of the detectors when they pass through a bright 
source. Our next step would be to try to understand this drifting behavior, and 
to incorporate it into our algorithm. 

CONCLUSION 

Our algorithm is just another way of trying to make use of the redundancy in the 
measurements. The standard HiRes processing uses it to pull out high spatial 
frequency information, while we take it a step further to let the scanlines cali- 
brate each other. It resides in the same theoretical framework as the Maximum 
Correlation Method. 
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